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This  report  pres— ts  the  results  of  —  AF0SR  spon— red  project  to  model  the 
homogeneous  gas  phe—  oxidation  kinetics  for  boron  — slsted  hydrocarbon  combustion  and  u— 
kinetic  s— sitivity  analysis  to  determine  critical  ruction  pathways  end  as— ss 
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In  addition,  a  preliminary  aval ua cion  of  cha  modifications  and/or  extensions  to 
sensitivity  analysis  techniques  for  future  applications  to  the  heterogeneous  phaaes 
of  boron  combustion  waa  performed. 

The  basic  foundation  has • been  developed  for  Initiating  experimental  and/or 
theoretical  studies  to  further  refine  the  kinetic  and  tharnochemical  parameters 
Identified  to  be  of  critical  Importance  to  gas  phase  oxidation  and  to  help  validate 
the  present  gas  phase  combustion  modal.  At  the  same  time,  a  fairly  comprehensive 
gas  phase  combustion  model  nov  exists  for  studying  a  number  of  boron  based  combustion 
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SUMMARY 


Solid  boron  ha*  long  boon  recognized  as  a  potential  high  volumetric 
density  propellant  for  volume  Halted  air  breathing  propulsion  systems* 
However,  achieving  this  potential  has  proved  difficult  due  to  low  combustion 
and  heat  release  efficiencies  encountered  when  solid  boron  slurries  are  added 
to  conventionally  design  coabustors*  The  difficulties  in  achieving  the 
anticipated  potential  of  boron- assisted  propellants  in  practical  applications 
has  created  a  need  for  laproved  understanding  of  the  basic  chealcal  and 
physical  processes  In  these  systems* 

The  coabustlon  of  solid  boron  particles  with  an  oxygen  based  oxldlser  has 
long  been  known  to  be  a  complex  process  involving  both  heterogeneous 
(gas/ solid)  and  homogeneous  (gas/gas)  chealcal  kinetics.  The  coabustlon 
process  begins  with  an  Initial  heterogeneous  phase  of  particle  oxidation. 
Ignition  and  volatilisation  of  partially  oxidised  boron  species.  The 
volatilised  boron  species  undergo  further  homogeneous  gas-phase  oxidation  to 
fora* a  mixture  of  gaseous  boron  oxide  and  boron  oxyhydride  species.  Finally, 
as  these  coabustlon  gases  cool, the  coaplex  B-O-H  mixture  must  condense  to 
liquid  or  solid  l20j  In  a  third  (heterogeneous)  coabustlon  phase  In  order  to 
coaplete  the  heat  release  necessary  for  an  efficient  boron  assisted 
propellant. 

This  study  presents  the  results  of  an  AFOSR  sponsored  project  to  aodel 
the  homogeneous  gas  phase  coabustlon  sub-process  and  use  kinetic  sensitivity 
analysis  to  determine  critical  reaction  pathways  and  uncertainties  In  the 
currently  available  theraochealcal  and  kinetic  data.  An  additional  objective 
waa  to  evaluate  the  extensions  of  sensitivity  analysis  that  would  permit  its 
future  application  to  heterogeneous  boron  coabustlon.  Eventually,  the  nodels 
for  each  of  the  three  sub-processes  would  be  coupled  to  develop  a  coaplete 
aodel  for  the  full  coabustlon  process. 
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From  the  modeling  work  and  sensitivity  analysis  summarized  in  this 
report,  several  important  results  hare  bean  obtained.  These  include  the 
followings 

1)  Thermodynamically,  the  dominant  equilibrium  boron  species  in  the 
boundary  layer  surrounding  a  hot  boron  particle  in  hydrocarbon 
combustion  gases  consists  of  HBO  and  B202  for  fuel  rich  systems  and 
U02  and  B203  for  oxygen  rich  systems. 

2)  The  initial  identity  of  boron  sub-oxides  and  oxyhydride  species 
mind  from  the  particulate  boron  Influences  the  gas  phase 
oxidation  kinetics  only  during  the  first  few  micro  seconds.  The 
species  snd  temperature  profiles,  as  well  as  sensitivity  gradients 
are  only  weakly  dependent  on  the  identity  of  initially  vaporized 
species  during  the  subsequent  oxidation  chemistry. 

3)  The  preaence  of  hydrogen  containing  species  has  a  significant  Impact 
on  the  homogeneous  combustion.  In  general  hydrogen  tends  to  have  a 
favorable  effect  because  it  aids  vaporisation  of  B(s)  and  promotes 
reactions  of  BO  snd  B02  to  fora  HB02  and  B^j.  On  the  other  hand, 
larger  quantities  of  hydrogen  promote  the  formation  of  HB02,  and  the 
product  condensation  stags  becomes  more  complicated. 

4)  The  gas  phase  oxidation  kinetics  is  dominated  by  a  relatively  small 
set  of  kinetic  and  equilibrium  processes.  These  have  been 
identified  for  further  experimental  and/or  theoretical  study. 

5)  HB02  and  B202  era  the  two  principal  gas  phase  combustion  products 
with  B203(g)  favored  at  higher  combustion  temperatures.  For  nost 
systems  modeled,  HB02  la  favored.  Consequently,  the  product 
condensation  stage  of  boron  combustion  is  probably  more  typically  a 
complex  heterogeneous  process.  However,  the  ability  to  shift  the 
product  distribution  to  B203  via  gptclMl  additives  would  allow 
condensation  to  proceed  as  a  simple  nucleation  which  could 
significantly  Increase  combustion  efficiency. 
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6)  A  preliminary  analysis  has  identified  key  issues  and  viable 

approachea  for  extending  the  current  model  and  sensitivity  analysis 
to  the  heterogeneous  burn  phase  of  boron  combustion. 

The  basic  foundation  has  now  been  developed  for  initiating  experimental 
and/or  theoretical  studies  to  further  refine  the  kinetic  and  thermochemical 
paraaMters  identified  to  be  of  critical  Importance  to  gas  phase  oxidation  and 
to  help  validate  the  present  gas  phase  combustion  model.  At  the  same  time,  a 
fairly  comprehensive  gas  phase  combustion  model  now  exists  for  studying  a 
number  of  boron  based  combustion  systems  of  important  to  the  Department  of 
Defense  community. 
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•tvilu  to  focus  effort  on  quantifying  only  those  nodal  parameters,  Including 
reaction  rets  constants  end  thamocheni cal  paraaatars,  which  control  affactlva 
boron  contention.  This  process  is,  of  course,  an  itarativa  procedure  with 
nodal  predictions  end  sensitivity  analyses  beconing  nore  accurate  as  input 
dots  for  important  physical  and  chemical  paraaatars  are  improved. 

▲  crucial  set  of  model  parameters  for  hydrocarbon  assisted  boron 
combustion  is  the  gas  phase  reaction  rate  constants  governing  the  oxidation  of 
suboxide  end  suboxyhydrlde  species  volatilised  during  the  combustion  of  solid 
particulate  boron.  Chemical  equilibrium  calculations  predict  that  solid  boron 
burning  in  the  presence  of  hydrocarbon  fuels  will  evolve  subetantlally  to 
gaseous  boron  oxyhydrlde  compounds  such  as  BOBO.  Since  the  chemical 
environment  in  the  boundary  layer  around  each  boron  particle  will  almost 
certainly  be  oxidiser  poor,  it  is  likely  that  a  I/O  ratio  of  less  than  two 
will  prevail  in  the  vaporised  species.  Depending  on  the  currently  uncertain 
details  of  gas/ surface  thermochemistry  sad  kinetics,  it  is  expected  thet  the 
volatilisation  of  species  such  as  BO,  8B0,  BOB,  B2O,  and  B^  as  the 
boron/ boron  oxide  surface  layer  of  the  particle  is  heated  by  radiative  and 
convective  heat  transfer  and  reacts  with  hydrocarbon  combustion  Intermediates 
and  products  including  H,  OB,  B20,  00,  and  002. 

Tbs  initial  modeling  effort,  presented  in  this  report,  concentrates  on 
the  homogeneous  gas  phase  portion  of  boron  combustion  for  several  related 

reasons.  First,  an  extensive  gas  phase  reaction  rate  set  for  boron  assisted 

eeubuetlon  was  developed  by  Kolb  and  co-workers  under  DABPA  and  AFOSl  funding 
la  the  late  IWO'i,1"1*1’**1"4  and  this  reaction  set  with  its  estimated 
rate  constant  parameters  can  serve  as  the  basis  for  initial  model 
development.  Second,  while  a  few  key  reactions  in  the  gee  phase  boron 
oxidation  scheme  have  been  measured  over  the  peat  several  years,  (generally 
confirming  the  rate  estimates  of  Kolb  and  co-workers)1'1  and  several 
experimental  groups  have  the  capability  to  measure  nore,  it  is  clear  chat  a 
kinetic  sensitivity  analysis  study  of  the  systen  will  help  specify  which  of 

the  many  unmeasured  rate  constants  should  be  the  focus  of  dedicated 

experimental  or  theoretical  efforts.  Finally,  kinetic  sensitivity  studies  of 
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combustion  systems  have  so  far  been  limited  to  homogeneous,  gas  phase  systems; 
thus,  the  gas  phase  portion  of  the  boron  combustion  process  is  most  amenable 
to  immediate  kinetic  sensitivity  analysis. 

An  important  secondary  task  of  our  initial  effort  is  to  map  out  the 
theoretical  extensions  necessary  to  extend  the  gradient  kinetic  sensitivity 
analysis  technique  to  the  heterogeneous  problems  of  boron  particle  ignition 
and  burnout  and  boron  oxyhydrlde/ boron  oxide  condensation  kinetics  necessary 
for  a  complete  boron  combustion  model. 

Section  2  of  this  report  will  present  elements  of  the  overall  model  for 
boron  assisted  combustion  of  hydrocarbon  fuels.  Section  3  presents  the 
homogenous  combustion  model  developed  to  date  and  the  results  of  initial  model 
runs  and'v«ensltlvlty  analysis.  Section  4  summarises  conclusions  and  presents 
recommendations  for  further  work. 
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2.  COMBUSTION  MODEL  DEVELOPMENT  AND  SENSITIVITY  ANALYSIS 


2.1  Modal  Overvlei r 

Particulate  boron  combustion  lu  conjunction  with  hydrocarbon  coabuatlon 
la  a  coaplax,  multiphase  process.2-1  Initial  analysis  lndlcatas  that  the 
process  can  be  viewed  as  Involving  at  laast  three  stages.  These  are:  1) 
boron  particle  surface  oxidation.  Ignition  and  volatilisation  of  suboxides  and 
suboxyhydrldes;  2)  a  largely  hoaogeneous,  gas  phase  oxidation  of  the  volatile 
suboxides  and  suboxyhydrlde  which  pass  through  the  oxldlser  poor  boundary 
layers  of  the  particle;  and  3)  condensation  of  the  gaseous  high  tenperature 
boron  oxyhydride  products  to  liquid  borla  (B2O3)  as  the  conbustlon  exhaust 
tenperature  and  pressure  approach  anblent  values. 


A  schematic  of  these  conbustlon  stages  Is  shown  In  Figure  2.1.  Each  of 
these  Is  discussed  In  nore  detail  In  subsequent  subsections. 


COMBUSTION  PHASE 


DESCRIPTION 


TYPICAL  SPECIES 


I 


II 


III 


BORON  PARTICLE  OXIDATION, 
IGNITION  AND  VOLATILIZATION 


HOMOGENEOUS  REACTIONS  OF  BORON 
SUBOXIDES  AND  GKYHYDRIDES 


CONDENSATION  OF  BORON  QXYHYDRIDES 
TO  B203(I) 


B(s) 
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HBO(g) 
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BO(g) 

C02(g) 
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(HOBO), 

B203<g) 

HxOyB,(l) 

B203(1) 

Figure  2.1.  Schematic  of  Boron  Particle  Conbustlon  Processes 
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2*2  Particle  Ignition  and  Volatilisation 


Tha  ignition  process  involves  tha  heating  of  an  initially  cold  boron 
particle  covered  with  an  oxide  layer.  A a  tha  particle  tenperature  rises  (due 
initially  to  the  hot  external  environment  from  the  surrounding  hydrocarbon 
combustion,  but  eventually  to  soma  self-heating),  two  competing  processes 
occur  Involving  the  removal  of  and  the  formation  of  the  oxide  layer.  The 
oxide  layer  is  removed  at  the  surface  through  vaporisation.  At  high 
temperatures,  this  process  can  Involve  chemical  transformations,  and 
therefore,  the  composition  of  the  surrounding  gases  can  significantly  alter 
the  rates  of  axldelsyer  removal. 2~ 2  For  example,  meter  vapor  has  been 
observed  to  promote  ignition. 

Soma  controversy  has  recently  arisen  regarding  the  mechanism  for  oxide 
formation.  Early  works2"3  postulated  that  oxygen  dissolves  in  the  oxide  layer 
and  diffuses  to  the  boron  -  boron  oxide  surface  where  it  reacts  to  form  boron 
oxide.  More  recently, 2-4  it  has  been  argued  that  a  more  likely  mechanism 
would  Involved  the  diffusion  of  boron  to  the  boron  oxide  -  gas  surface  where 
subsequent  heterogeneous  reaction  with  oxygen  to  form  the  oxide  would  occur. 

Most  of  the  past  work  on  boron  particle  ignition  has  Involved  studies  of 
boron  in  oxygen  or  air.  However,  as  noted  above,  water  vapor  can  have  a 
significant  promotional  effect  on  boron  ignition  kinetics.2-2  In  a  boron 
assisted  hydrocarbon  combustion  process,  not  only  water  vapor  but  a  number  of 
other  hydrocarbon  combustion  intermediates  and  products  including  H2,  00,  OH, 
B,  C02,  and  possibly  partially  oxidised  hydrocarbons  such  as  aldehydes  will 
all  Interact  with  the  surface  of  the  boron  particles. 

It  is  probable  that  it  is  the  chemical  reaction  of  these  species,  not  0 
and  02  which  control  particle  ignition  and  vaporisation  in  mixed 
hydrocarbon/boron  combustion  systems.  Since  both  thermochemical 
models,2-5-2-6  and  flame  sampling  studies2-7  indicate  that  boron  oxyhydrlde 
species  such  as  HOBO  are  the  dominant  high  temperature  boron  containing 
combustion  products,  it  seems  possible  that  a  third  process  dominates  boron 
particle  ignition  and  volatilisation  kinetics  in  hydrocarbon  combustion 


iTTvv7yv\ 


ays t am .  This  procsss  would  Involve  Che  removal  of  Che  boron  oxide  layer 
itself  through  surface  reaction  with  hydrogen  containing  species  such  as  H, 

OH,  Hjt  and  H2O  to  fore  volatile  boron  oxyhydride  species  such  as  HBO,  BOH, 
and  HOBO.  As  the  depth  of  the  boron  oxide  layer  is  reduced,  new  surface  will 
be  available  to  refora  boron  oxides,  through  interaction  with  002,  OH,  H20, 
and  whatever  0  and  02  penetrate  the  oxygen  poor  region  around  the  particle. 

Indeed  the  reactive  chemistry  of  atomic  boron  as  well  as  electronically 
analogous  atomic  aluminum  provides  a  strong  hint  that  oxyhydride  species  form 
efficiently  and  preferentially.  These  electron  deficient  atoms  each  combine 
with  water  in  bimolecular  reactions  in  both  the  gas  phase2-8-2"*12  and  in  cold 
argon  matrices2-13  to  form  an  excited  H-M-OH  species  which  efficiently  sheds 
its  excess  energy  through  radiative  decay. 

B,  A1  +  H20  ♦  H-B-OH*  or  H-A1-0H*  -*•  H-B-OH  or  H-Ai-OH  +  hv  (2.3.1) 

Theoretical  studies  confirm  that  these  remarkable  reactions  occur  with 
surprisingly  small  activation  energies.2-14-2-15  It  is  not  difficult  to 
imagine  that  the  high  combustion  temperatures  coupled  with  the  presence  of 
surface  to  stabilise  reactions  may  well  promote  very  efficient  oxyhydride 
formation. 

As  surface  reactions  and  heat  transfer  raise  the  particle  temperature, 
boron  suboxide  species  such  as  BO,  B20,  and  B202  may  vaporize  along  with  the 
boron  oxyhydride  species. 

Thus,  it  is  quite  possible  that  under  hydrocarbon  combustion  conditions 
the  boron  oxide  layer  which  strongly  inhibits  particle  ignition  and 
vaporisation  under  clean  boron/oxygen  system  conditions  becomes  a  reactive 
layer  promoting  particle  vaporization  when  the  particle  boundary  layer  is 
oxidizer  poor. 


A  tatter  understanding  of  the  thereochenical,  solubility ,  transport,  and 
reaction  properties  of  the  layer  is  required  to  resolve  the  actual  ignition 
end  vaporisation  aechanlsu.  Thermodynamically,  the  oxide  formation  process  is 
exothermic  whereas  the  removal  process  for  oxides  is  endothermic.  Thus,  the 
criteria  for  Ignition  in  the  "clean"  boron/oxygen  system  Is  that  the  oxide 
vaporisation  rate  exceeds  its  production  rate  so  that  the  layer  Itself  Is 
removed.  The  degree  to  which  hydrocarbon  combustion  products  and 
intermediates  change  this  picture  is  an  Important  question. 


itlon  Kinetics 


Once  the  oxide  layer  is  removed,  the  particle  burns  rapidly  with  a 
surrounding  flame  sheet  which  is  detached  from  but  near  the  surface.  Since 
the  volatility  of  boron  is  low,  boron  burns  much  like  a  carbon  particle  vs.  a 
hydrocarbon  liquid  droplet. Here,  gasification  of  the  boron  occurs  by 
chemical  transformation  into  more  volatile  components.  Thus,  heterogeneous 
end  homogeneous  chemistry  play  a  more  dominant  role  in  the  gasification  of 
boron  than  for  liquid  hydrocarbons.  For  example,  liquid  hydrocarbon  droplets 
gasify  primarily  by  a  simple  vaporisation  process  which  proceeds  in  the 
presence  or  absence  of  oxidant.  The  volatile  component  formed  during  boron 
gasification  is  believed  to  be  BO,  but  as  noted  above  there  may  be  significant 
levels  of  HOBO,  HBO,  BOH,  B2O,  or  BjOj.  C02  is  the  major  product  of  carbon 
combustion  which  diffuses  back  to  the  carbon  particle  surface  to  react 
heterogeoeously  to  form  more  CO.  Boron  oxides  and  oxyhydrldes  may  play 
similar  roles  in  boron  coabustlon.  Indeed,  there  nay  be  a  strong  similarity 
between  boron  and  carbon  combustion;  even  the- homogeneous  gas  phase  kinetics 
may  have  similarities.  For  example,  CO  oxidation  is  controlled  by  the 
presence  of  the  hydrogen  containing  OH  radical  and  boron  sub-oxide  oxidation 
will  be  clearly  affected  by  oxyhydrlde  formation. 


Evolution  of  gas-phase  boron  compounds  in  the  boundary  layer  around  the 
burning  particle  is  most  likely  under  kinetic  control  and  must  be  modeled  as  a 
homogeneous  kinetic  process.  This  requires  estimated  or  measured  rate 
constants  for  a  large  number  of  individual  reactions.  Fortunately,  previous 


DOD  interest  In  boron  assisted  liquid  rocket  propellants  and  their  exhaust 
pluses  resulted  in  considerable  activity  in  this  area.2-16  Section  3  presents 
an  updated  and  extensive  version  of  this  reaction  rate  constant  set, 
including,  of  course,  the  results  of  available  experimental  data  obtained 
since  1979. 

2.4  Product  Condensation 

Estimates  for  the  high  volumetric  energy  density  potential  of  boron-based 
fuels  are  based  on  heats  of  combustion  for  complete  oxidation  with  condensed 
boron  oxide,  e.g.,  B203  (1),  as  the  principal  product.  Thus,  realising  the 
full  theoretical  potential  depends,  in  part,  on  achieving  complete  energy 
conversion  as  combustion  product  vapors  cool  during  the  expansion  process. 

This  has  been  difficult  to  achieve,  however,  because  chemical  recombination 
reaction  ratea  are  typically  slow  relative  to  the  short  residence-times  for 
expansions  in  typical  combustion  devices.  This  significantly  reduces 
combustion  efficiency  as  energy  is  effectively  trapped  in  the  vapor  phase.  It 
is  currently  estimated  that  as  much  as  20  -  30  percent  of  the  estimated 
available  energy  nay  be  lost  by  this  mechanism.  The  problem  is  particularly 
acute  for  high  altitude  operation  where  reaction  rates  are  even  slower  due  to 
the  reduced  pressure  in  the  nossle  exhaust  flow. 

Equilibrium  calculations  for  boron  oxidation  in  wet  flames  Indicate  large 
concentrations  of  HOBO  in  the  product  vapors.  Studies  of  the  boron-oxide 
surface  layers  properties  also  suggest  the  heterogeneous  formation  of  HOBO 
during  combustion  in  the  presence  of  H20.  Taking  H0B0(g)  as  the  principal 
vapor  phase  product,  boron  oxidation  at  temperatures  near  the  onset  of 
nudeation  (T  '2000  K)  may  be  represented  as 


2B  +  02  +  H2  ♦  2H0BO  (g) 


AH  -  1151  KJ/mol 


(2.4.1) 


In  contrast,  complete  combustion  to  condensed  boron  oxide. 


2B  +  202  +  2  ♦  H2  ♦  B^U)  +  H2°  48  "  1468  K^*01  (2.4.2) 

is  approximately  20  percent  more  exothermic.  Since,  at  T  *  2000  K,  the 
constant  pressure  equilibrium  constants  for  reactions  2.4.1  and  2.4.2  are  to  a 
first  approximation  equal,  failure  of  H0B0(g)  to  convert  to  B202  and 
subsequently  condense  significantly  reduces  combustion  efficiency.  Equations 
(2.4.1)  and  (2.4.2)  do  not  represent  the  reaction  mechanisms  for  boron 
oxidation,  but  only  serve  to  indicate  the  Increased  energy  available  upon 
condensation.  Obtaining  this  extra  energy,  however,  is  dependent  upon 
condensation. 

The  mechanism  for  formation  of  B203(1)  remains  one  of  the  most  poorly 
understood  components  of  boron  combustion.  Edelmann  et.  al.  characterized 
it  as  a  nudeatlon  process  and  developed  a  model  for  the  finite-rate  kinetics 
and  condensation  during  nossle  expansions  which  agreed  reasonably  well  with 
experiment.2-17  In  contrast.  Miller  hss  argued  that  the  mechanism  is  more 
complex  than  just  a  simple  nudeatlon  process  and  actually  Involves  many 
intermediates  and  reactions  not  considered  in  Edelmann' s  model.2-7  Miller's 
arguments  are  based  on  studies  of  the  boron  combustion  products  at  low 
pressures,  ca.  0.1  atm,  in  which  mass  spectrometry  was  used  to  determine  the 
species  present  in  sn  oxygen-rich  fleam  under  conditions  appropriate  for 
nudeatlon.  These  results  indicated  the  presence  of  a  variety  of  complex 
boron  oxyhydroxlde  species  in  addition  to  simple  monomers,  dimers  and  trimers 
of  HOBO.2-7  Since  polymerisation  reactions  are  typically  fast  (k  ”10_1°  cm3 
molecule  sec  for  dimer  and  trlmer  formation),  Miller  concluded  that 
elimination  of  H20  from  the  polymer  chain  was  the  more  likely  kinetic 
bottleneck  to  condensation.  Moreover,  it  appeared  that  H20  elimination 
occurred  gradually  during  nudeatlon  rather  than  in  a  single  discrete  step. 


From  this  data.  Miller  formulated  a  hypothetical  polymerization 
mechanism.  However,  the  experimental  studies  were  unable  to  probe  specific 
reaction  pathways,  and  many  reasonable  pathways  were  left  out  of  the  model,  so 
that  no  quantitative  predictions  that  could  be  compared  with  experiment  were 
made.  It  has  been  suggested  that  some  of  the  discrepancy  between  Edelman  et. 
al.  and  Miller  may  be  attributed  to  the  fact  that  they  had  examined  the 
combustion  products  in  very  different  pressure  ranges.  However,  since  neither 
study  provided  definite  measurements  to  test  reaction  mechanisms,  both  models 
are  subject  to  a  great  deal  of  uncertainty.  Since  the  formation  of  condensed 
boron  oxide  Is  so  Important  to  combustion  efficiency,  it  Is  imperative  to  gain 
a  more  reliable  grasp  of  the  key  reactions  and  parameters  which  control  this 
process. 
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3.  HOMOGENEOUS  BORON  COMBUSTION:  MODEL  DEVELOPMENT  AND  KINETIC 
SENSITIVITY  ANALYSIS 

3.1  Overview  of  Model  Development 

The  model  development  and  concomitant  sensitivity  analysis  program 
described  in  the  report  consisted  of  the  following  tasks: 

1)  The  collection  snd  evaluation  of  currently  available  thermochemical 
and  kinetic  data  for  chemical  species  and  elementary  reactions  that 
may  be  Important  to  boron  combustion  (Subsection  3.2). 

2)  A  series  of  thermodynamic  calculations  for  air  oxidation  of  a 
mixture  of  solid  boron  and  a  typical  hydrocarbon  fuel  (JP4)  to 
determine  equilibrium  boron  speciation  (Subsection  3.3). 

3)  The  development  of  a  homogeneous  combustion  model  consisting  of  a 
detailed  list  of  plausible  elementary  reactions  and  estimated  (or 
measured  values  when  available)  rate  constants  followed  by  kinetic 
calculations  to  determine  species  concentrations  and  temperature 
profiles  under  various  Initial  conditions.  (Subsection  3.4). 

4)  The  application  of  kinetic  sensitivity  analysis  to  determine  the 
response  of  system  observables  to  changes  in  the  input  parameters, 
e.g..  Initial  species  concentrating,  rate  parameters,  etc. 

5)  An  In  depth  analysis  of  the  gradient  sensitivity  data  to  identify 
the  critical  reaction  pathways  and  key  uncertainties  in  currently 
available  thermochemical  and  kinetic  parameters  (Subsection  3.6). 

6)  and  lastly,  a  brief  analysis  of  the  research  strategies  needed  to 
extend  the  applicability  of  sensitivity  analysis  to  heterogeneous 
process. 

Each  phase  of  this  work  Is  briefly  described  In  the  following  subsections. 
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3.2  Collection  and  Evaluation  of  talevant  Thai 


Kinetic  Dote 


rmnrhamlcal  and 

During  the  initial  stage  of  this  program,  a  eoaprehenslvt  Literature 
search  was  performed  to  develop  a  data  base  containing  tha  most  up-to-date 
values  for  the  thermochemical  and  kinetics  parameters  that  would  he  needed  In 
developing  a  nod el  for  gas  phase  boron  oxidation  kinetics.  The  data  bene  woe 
periodically  re-evaluated  and  updated  as  new  or  acre  reliable  results  became 
available.  In  this  section,  we  suamarlse  the  more  Important  results  from  this 
survey. 


3.2.1  Thermochemical  Data 

Table  3-1  contains  the  standard  enthalpies  of  forantiom  &df,2M 
entropies  S298  for  the  atonic  and  molecular  species  that  are  lnportaat  to 
boron  assisted  hydrocarbon  combustion  ays tame .  (This  list  intentionally 
concentrates  on  species  containing  boron  since  the  corresponding 
carbon-containing  species  are  well-known  and  well-doc  wanted . )  With  one 
notable  exception  (BOH),  this  data  has  been  collected  from  the  meet  recent 
publication  of  the  JANAF  thermochemical  tab las. 3 The  heavier  boron 
oxyhydrldes  (e.g.  H3B03,  H^BjO,,,  etc.)  are  probably  not  critical  to  a 
description  of  the  high  temperature  gee  phase  boron  coobmetlom  chemistry. 

They  will  be  important,  however,  to  future  modeling  of  the  condensation  phase 
(Subsection  2.3)  and  are  included  in  Table  3-1  for  completeness. 

The  HBO  species  has  long  been  thought  to  be  an  lnportaat  Intermediate  In 
boron  conbustlon  system  containing  hydrogen.  It  Is  pertinent,  therefore,  to 
consider  the  possible  role  of  the  BOH  Isomer.  At  present  there  Is  no 
experimental  data  available  for  this  species.  Ab  initio  calculations, 
however,  have  been  performed  to  characterise  the  potential  energy  and  mini mum 
energy  path  for  the  intramolecular  rearrangement  BOH  *  IBO.3-2  These  results 
Indicate  that  the  ground  configuration  for  BOH  is  approxlnstely  40  kcal  above 
the  ground  configuration  for  HBO.  The  potential  barrier  for  rearrangement  is 
approximately  SO  kcal.  The  optlnun  geometric  perameters  for  the  BOH  species 
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Tram  Tabla  3-1 ,  it  la  lanediately  apparent  chat  cba  uncertainties  In  the 
tharaMpuMe  parameters  for  tha  B/O/H  combustion  ays  tea  are  several  orders  of 
nagnitads  larger  than  for  tha  sane  paranatars  associated  with  nolacular 
apaalaa  lapsrtaat  to  conventional  hydrocarbon  combustion.  For  example »  the 
ancartalntlaa  In  the  standard  enthalpies  and  entropies  for  CO  and  C02  range 
from  ana  to  tmo  orders  of  magnitude  smaller  than  tha  corresponding 
umeertalntles  their  laoalactronle  boron  analogs  HBO  and  HB02 . 

The  heats  of  fonsstlon  In  Tabla  3-1  have  bean  used  to  establish  bond 
energies  far  several  relevant  1-H,  B-0,  and  B-B  bonds.  These  arc  sunnarlsed 


la  Table  3-2. 


Table  3-2  -  Bond  Inarglea  for  Boron  Assisted  Hydrocarbon  Combustion 
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118  ±  0.1 
256  ±  0.04 

127.6  ±  0.06 
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3 .2 .2  Experimental  Kinetic  Data 


Reliable  data  for  kinetic  parameters  associated  with  the  B/O/H  combustion 
system  is  even  more  lacking  than  for  the  thermodynamic  parameters.  One 

important  exception  is  the  recent  work  by  Baughcum  and  Oldenborg  at  Los 

Alamos  National  Laboratory  on  the  reactions  of  B  and  BO  with  molecular 
oxygen. 

B  +  02  -►  BO  +  0  kf  -  1.19  (±0.004)  x  1010  e-(158  *  13)/T  cm3/sec 

BO  +  02  B02  +  0  kf  -  1.7  (±0.4)  x  10_1  *(1/1000)  2.23  ±  0.5 

(716  ±  288)/T  3. 

e-  cm  /sec 

This  work  extends  the  earlier  room  temperature  work  on  B  +  02  by  DIGiusspe  and 
Davidovits3-**  and  BO  +  O2  by  Llewellyn,  Fontijn  and  Clyne.3-5  We  have  also 
drawn  on  the  experimental  work  of  DiGiusseppe  et.  al.,3-i*  "  3-6  for  the 
reaction  of  B  with  002  and  B  with  H20. 

3.3  Combustion  Thermochemistry 

Equilibrium  calculations  were  performed  to  establish  Initial  estimates  of 
the  boron  speciation  in  the  boundary  layer  surrounding  each  boron  particle. 

The  results  of  these  calculations  were  used  to  extend  the  list  of  relevant 
reaction  (Subsection  3.4.1)  and  as  guides  in  choosing  initial  species 
concentrations  for  the  subsequent  kinetic  calculations  (Subsection  3.4.3). 

The  equilibrium  calculations  were  performed  using  the  PACKAGE  code 
(Plasma  Analysis,  Chemical  Kinetics,  and  Generator  Efficiency)  developed  at 
ARI  to  predict  multiphase  combustion  product  parameters  including  mixture 
composition  and  thermodynamic  properties.  Originally  based  on  the  NASA  CED 
code,  the  present  PACKAGE  program  has  been  extended  to  Include  a  multiple 
ideal  solution  model,  a  phase  rule  check  when  condensed  phases  are  present, 
and  improved  methods  for  dealing  with  numerical  difficulties. 
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The  equilibrium  calculations  wars  parforaad  for  oxidation  of  a  JP4/ solid 
boron  aixtura  (foal)  in  air  (oxidisar).  Initial  concantrationa  of  tha  fual 
roaetanta  (aolaa  of  fual  raaetant/ total  aolas  of  fual)  and  of  tha  components 
of  tha  oxidising  medium  (aolas  of  oxidisar  spaciaa/total  aolaa  of  oxidisar) 
ara  givan  in  Tabla  3-3  and  wars  tha  aaaa  for  aaeh  calculation.  The  actual 
raaetant  concantrationa  in  tha  fual-oxldlsar  aixtura  vara  varied  by  adjusting 
tha  fual/ oxidisar  (F/0)  equivalence  ratio.  Calculationa  wars  parforaad  for 
7/0  equivalence  ratios  ranging  froa  0.3  (oxidisar  rich)  to  4.3  (fual  rich). 
However,  since  tha  chemical  environment  in  tha  boundary  layer  is  almost 
certainly  fual  rich,  the  calculations  for  7/0  >2.0  are  probably  aore 
representative  of  the  boron  species  concantrationa  near  a  boron  particle. 

Tabla  5-3  -  Initial  Tual  and  Oxidisar  Species  Concentrations  for 
Iquillbritss  Calculationa 


Species 

Hole  Fraction 

Fuel: 

JP4 

0.1 

»<e> 

0.9 

Oxidisar : 

*2 

0.78 

°2 

0.21 

Ar 

0.0095 

co* 

0.0005 

Figures  3.1a  and  3.1b  show  constant  pressure  equilibrium  species 
concentrations  as  a  function  of  temperature  for  7/0  equivalence  ratios  of  0.3 
and  4.3,  respectively.  Tha  pressure  was  1.0  atn  end  the  tenperature  was 
varied  froa  1400  K  to  2400  K.  Similar  results  for  pressures  of  3.0  atm  and 
6.0  atn  did  not  indicate  any  substantial  change  in  the  thermodynamically 
predicted  boron  sped  at  ion. 
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Figure  3.1b.  Equilibria*  Concentration  Versus  Teaperature  for  Boron 

Spaciatlon  in  tha  Boundary  Layer  of  a  Hot  Boron  Particle 
in  Hydrocarbon  Coabustlon  Gases.  (P  -  1.0  ata;  F/O 
Equivalence  Ratio  *4.3  ata) 
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For  F/0  •  0*3,  Figure  3.1a,  the  system  is  oxidizer  rich  end  the  Initial 
aolld  boron  phaae,  excluding  the  liquid  borla  component  [B203(1)]  below  1800 
K,  la  completely  oxidized  to  gas  phase  oxides  and  oxyhydrides.  The  dominant 
gas  phase  species  are  HOBO,  B203,  and  B02  for  T  >  1600  K.  The  concentrations 
for  the  heavier  oxyhydrides  H3B03  and  H3B306,  although  comparable  to  the 
concentrations  of  B203  and  B02  for  T  <  1600  K,  are  more  than  3  orders  of 
magnitude  smaller  then  the  concentration  of  HOBO  for  T  “  1800  K.  Thus, 
although  these  species  will  be  important  in  modeling  the  condensation  stage  of 
boron  combustion,  they  are  not  going  to  be  important  to  the  gas  phase  kinetics 
considered  in  this  report.  Beginning  at  about  2200  K,  entropy  terms  begin  to 
become  important  and  there  is  a  gradual  increase  in  the  concentrations  of 
smaller  species,  B0  and  HB0,  with  Increasing  temperature.  However,  they 
remain  minor  species  over  the  temperature  range  considered  with  concentrations 
at  least  3  orders  of  magnitude  smaller  then  the  the  dominant  species  HOBO. 

The  equilibrium  boron  epecies  distribution  changes  markedly  with  as  the 
F/O  equivalence  ratio  increases.  Figure  3.1b  shows  the  species  concentrations 
for  an  equivalence  ratio  of  4.5.  These  results  are  typical  for  fuel  rich 
systems  (F/O  >  2.0).  In  this  case,  a  solid  boron  phase  persists  throughout 
the  temperature  range  considered.  In  contrast  to  oxidizer  rich  systems,  the 
primary  equilibrium  gas  phase  species  for  fuel  rich  systems  are  HBO  and  B202. 
Gaseous  B203  and  HOBO  are  comparatively  minor  species  with  concentrations 
approximately  2  orders  of  magnitude  smaller. 

Table  3-4  summarizes  the  dominant  equilibrium  boron  speclation  in  the 
boundary  layer  surrounding  a  hot  boron  particle.  These  results  show  an 
interesting  correlation  with  conventional  hydrocarbon  combustion. 

Specifically,  the  dominant  combustion  products  for  boron  combustion,  HB02 
(oxidizer  rich)  and  HBO  (fuel  rich),  are  lsoelectronlc  with  the  primary 
hydrocarbon  oxidation  products  C02  (complete  combustion)  and  CO  (incomplete 
combustion),  respectively.  Similarly,  the  secondary  equilibrium  constituents 
show  expected  trends  in  the  B/0  mole  ratio  as  the  F/0  equivalence  ratio 
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Table  3-4  -  Dominant  Equilibrium  Gas  Phase  Boron  Speclation  at  T-1800  K 


P/0  Equivalence  Ratio: 

0.3 

Species 

Mole  Fraction 

Species 

Mole  Fraction 

HOBO 

10”2 

HB0 

10-2 

B203 

10“3 

B  2®  2 

10“  3 

bo2 

10“4 

B2O3 

10“4 

H0B0 

10“4 

increases.  Thus,  for  oxidizer  rich  systems,  the  B2O3  species  has  a  B/0  mole 
ratio  of  2/3,  whereas  for  fuel  rich  systems,  the  second  major  constituent  B202 
has  a  B/0  mole  ratio  of  1 . 

3.4  Homogeneous  Combustion  Kinetics 

3.4.1  Elementary  Reactions  and  Rate  Parameters 

There  is  currently  not  enough  known  about  the  initial  ignition  and 
volatilization  phase  of  boron  assisted  combustion  to  accurately  characterize 
the  boron  species  distributions  dominating  the  homogeneous  kinetics.  Kinetic 
and  transport  effects  may  significantly  alter  the  predicted  species 
concentrations  based  on  an  assumption  of  local  equilibrium.  The  equilibrium 
results  do  serve  as  guides,  however,  which  help  identify  the  species  and 
elementary  reactions  which  should  be  included  in  a  homogeneous  boron 
combustion  model.  Thus,  from  the  list  of  boron  species  in  Table  3-1,  the 
small  subset  given  in  Table  3-5  has  been  selected  for  initial  study.  The 
principal  species  not  included  are  the  heavier  boron  oxyhydrides  (HxBOy 
with  x,y  >  2)  which  appear  to  only  be  important  for  temperatures  below  1600  K 
and  those  species  (B20  and  BOH)  which  because  of  their  large  positive 
enthalpies  of  formation  tend  to  be  important  only  in  extremely  high 
temperature  regimes,  i.e.  T  >  2400  K. 

Table  3-5  -  Boron  Species  List  for  Homogeneous  Combustion  Models 
t - 1 


Based  on  Che  thermodynamically  predicted  boron  speciation  in  the  boundary 
layer  a  list  of  elementary  reactions  for  boron  assisted  hydrocarbon  combustion 
was  developed.  This  reaction  list,  which  is  summarized  in  Table  3-6,  includes 
most  of  the  thermodynamically  allowed  reactions  between  boron  containing 
species  and  hydrocarbon  intermediates.  Reactions  which  were  spin  forbidden, 
sterically  improbable,  or  not  thought  to  be  fundamental  were  not  included  in 
Table  3-6. 

Developing  an  accurate  model  for  the  homogeneous  boron  oxidation  kinetics 
requires  a  knowledge  of  the  rate  parameters  for  a  large  number  of  elementary 
reactions.  Unfortunately,  with  but  a  few  important  exceptions,  detailed 
experimental  and/or  theoretical  studies  are  notably  lacking  for  the  reactions 
listed  in  Table  3.4.  This  has  forced  this  initial  modeling  effort  to  rely  on 
estimated  rate  constants. 

Rate  constants  for  the  bimolecular  rate  constants  were  estimated  using 
transition  state  theory  (TST)  and/or  by  scaling  the  rate  constants  for  similar 
reactions.  Both  approaches  are  subject  to  large  uncertainties.  For  TST, 
general  questions  concerning  the  validity  of  a  statistical  theory  are 
relatively  minor  for  the  present  application;  rather,  the  primary  source  of 
uncertainty  stems  from  not  knowing  the  relevant  transition  state 
configurations  and  activation  energies.  Rate  constants  based  on  scaling  the 
rate  parameters  for  similar  reactions  do  not  necessarily  fare  much  better 
since  the  specific  scaling  relationship  used  may  also  require  either 
structural  and  dynamical  information  which  is  unknown,  or  be  forced  to  rely  on 
unsupported  assumptions.  Because  of  these  inherent  uncertainties,  the 
principal  aim  has  been  to  establish  rate  constants  which  are  at  least 
physically  reasonable  and,  to  a  certain  extent,  collectively  self-consistent. 

For  the  sets  of  reactions 

A  +  HBO  +  AH  +  BO  A-(H,  0,  and  OH)  (3.1) 

and 

A  +  HBO  H  +  ABO  A  -  (0,  OH,  and  BO)  (3.2) 
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HBO  +  OH- HO BO  +  H  -1.94  8x10 


the  rate  constants  wars  first  eetlaated  using  TtT  by  aeauniag  a  linear 
transition  stats.  Tbs  rats  constants  for  tbs  reactions  (3.2)  vara 
subsequently  decreased  by  about  an  order  of  aagnitnde  to  offset  tbe  affect  of 
assuaing  a  linear  coupler.  Kate  constants  for  reactions  (3.1)  vara  also 
calculated  by  scaling  the  corresponding  reactions  for  the  leoelectronic  BCM 
using  the  ratio  of  classical  partition  functions  for  the  reactants.  These 
values  were  typically  larger  than  the  rate  constants  calculated  aalag  TST. 
Since  TST  tends  to  yield  rate  constants  which  are  slower  than  experimentally 
observed,  the  TST  rate  constants  for  both  sets  of  reactions  ware  increased 
8 lightly  to  give  the  final  values  in  Table  3.4. 

Rate  constants  for  the  reactions 
B202  +  0  ♦  BO  +  IOj 

and 

B2 02  +  OH  -*•  BO  +  BOBO 


were  more  difficult  to  estinate  because  of  they  Involve  5  and  6  aton 
transition  states.  Initial  estinates  ware  obtained  by  slnply  doubling  the 
rate  constants  for  the  corresponding  HBO  reactions  (3.2)  due  to  the  syanetry 
of  the  B202  species  giving  kf  -  1.6xl0"u.  This  value  was  then  scaled  by 
the  reactant  partition  function  ratio  for  the  HBO  and  B202  reaction  giving  the 
value  in  Table  3-4,  kf  -  6x1 0"1 1 .  The  reactions 


HBO  +  OE  -*•  B02  +  Hj  kf 


-20  -35200/T 

10  xe 


HBO  +  02  +  B02  +  OH  kf 


,n-20  -35200/T 

10  xa 


B202  +  OH  ■*>  B02  +  HBO  kf 


-  10 


-15 


-35200/T 

xe 


B2°2  +  °2 


BO_  + 


B°2  kf 


,a-15  -35200/T 

■10  xe 


probably  involve  4-center  reaction  coaplexes  and,  consequently,  were  eetlaated 
to  proceed  very  slowly. 
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slgktai  spacific  baacs  ac  caastaae  precede  aai  tbs  aaaa  eelghcei 
for  aaab  species  are  gft—  by. 
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The  nIwmm  Mali  (pacific  haats  at  constant  pressure,  enthalpies,  and 
satrapies  far  each  species  ara  obtalaad  froa  tha  NASA  polynoainal  fits  to  the 

JUtf  data. 
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Tha  ala—  tary  reaction  rata  constant  la  expressed  la  aodlfiad  Arrhenius 


farm. 
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where  Afc  la  the  pre-expoanntlal  factor,  the  teaperstura  exponent,  and 
Ig'lt  tha  activation  energy.  f\  la  tha  cone t ant  parameter  perturbed  In  the 
seaeltivlty  aaalyala,  ^  for  each  reaction  equals  unity.  The  backward  rate 


constant  la  obtained  froa  tha  forward  reaction  rate  co net ant  and  the  reaction 
equlllbrlun  constant  using  tha  relation 


(3.4.15) 
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The  equlllbrlua  cone tents  for  the  reactions  ere  obtained  from 
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where  Kp,k  is  determined  froa 
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where 
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All  teraoleculer  reactions  ere  assuasd  to  be  in  their  low  pressure 
Halting  reglaes.  The  third  body  Mfe  for  each  reaction  is  defined  as  the  sua 
of  the  product  of  all  species  concentrations  tlaes  the  corresponding  chaperon 
efficiency 
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where  e^,!  is  the  chaperon  efficiency  for  the  1-th  species  and  the  k-th 
reaction  end  ci  is  the  concentration  of  the  i-th  chealcal  species. 
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The  governing  (quatlou  for  Cho  isothermal,  constant  volume  system  are 
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The  species  concentrations  are  related  to  tha  spades  mass  fractions  and  mole 
fractions  through  the  following  equations 


*  T  (  l  ) 


PWY. 

S"" pYi 


i-1  wi 


xi  • 


*  T. 

W4  (  I  ^  ) 
1  i-1  "i 


(3.4.22) 


(3.4.23) 


Tha  system  equations  are  solved  using  L80DK  and  several  modified 
subroutines  from  CUIKI*. 


3.4.3  jesfgnaf  of  8 secies  Concentration  and  Tan i^a return  Profile 

to  System  Far ama tars  ~ 

General  characteristic  behavior  of  high  temperature  l-C-H-0  chemistry  was 
investigated  by  examining  the  rosponso  of,  the  species  concentrations  and 
temperature  profiles  to  various  different  initial  conditions.  The  initial 
conditions  studied  covered  a  wide  range  of  temperature,  pressure,  and  initial 
species  concent rat ions.  These  conditions  wars  chosen  to  be  representative  of 
soma  of  the  gas  phase  environments  surrounding  metallic  boron  particle 
combustion.  Several  of  these  systems  are  Illustrated  below.  The  initial 
conditions  for  each  system  are  tabulated  in  Table  3-7. 

In  the  first  example  (tun  1),  the  initial  conditions  were  derived  from 
equilibrium  calculations  for  a  fuel  rich  mixture  of  JP4  and  B(s)  in  air  at  a  a 
pressure  of  1  atm  and  an  equivalence  ratio  of  4.5.  Results  from  these 


Tiklt  3*7  -  Initial  Coedit lees  for  tine tic  Calculations 
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eoaataat  pressure,  constant  temperature  cslculstions  over  a  tsapsraturs  range 
I  of  1400  -  2400  K  *re  shown  in  Figure  3.1b*  The  initial  conditions  for  the 

kinetic  calculations  warn  taken  at  a  t separatum  of  1800  K.  Here,  all  species 

■with  nole  frnctlons  greater  than  1  *  10“5  sera  considered  as  Initial 

reactants.  The  solid  boron  uns  replaced  by  an  equivalent  aaount  of  air  (where 
—  air  la  as sowed  as  212  O2  and  792  H2).  The  msults  of  species  nole  fractions 

|  and  tawperature  versus  tine  for  an  adiabatic,  constant  pressure  system  are 

shown  la  Figure  3.2.  Fart  (a)  of  the  figure  presents  the  boron  containing 
P  coup Quads;  part  (b)  presents  the  species  of  the  C0/H2/02  system;  and  part  (c) 

presents  the  tawperature. 
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Run  2,  Figure  3.3,  is  Che  sane  as  run  1  except  Chat  Che  metallic  boron 
removed  from  Che  system  in  run  1  was  not  replaced  with  any  species,  and  the 
existing  equilibrium  mole  fractions  were  normalized  to  equal  1.0,  neglecting 
solid  B. 

Runs  3  through  5,  Figures  3.4  through  3.6,  make  a  set  for  comparison. 

Run  3  is  for  an  initial  mixture  consisting  of  HBO,  H2,  CO,  02  and  N2  as 
reactants.  This  case  represents  an  extreme  assuming  all  the  initial  boron 
(e.g.,  coming  from  the  solid  boron  surface)  is  in  the  form  of  HBO.  Run  4  is 
the  same  as  run  3  except  the  reaction  is  Isothermal,  constant  pressure.  Run  5 
is  the  same  as  run  3  except  the  initial  H2  is  replaced  by  H20  and  the  initial 
CO  is  replaced  by  C02. 

Run  6-8,  Figures  3.7  through  3.9,  again  make  a  set.  Run  6  is  for  an 
initial  mixture  consisting  of  BO,  H2 ,  CO,  02  and  N2.  This  case  represents 
another  extreme  assuming  all  the  initial  boron  is  in  the  form  of  BO.  In 
comparison  to  run  3,  the  amount  of  BO  is  equal  to  the  amount  of  HBO.  The  H 
originally  in  HBO  (in  run  3)  is  put  into  H2  in  run  6.  Run  7  is  the  same  as 
run  6  except  the  reaction  is  Isothermal,  constant  pressure.  Run  8  is  the  same 
as  run  6  except  the  H2  is  replaced  by  H20  and  CO  by  C02. 

Run  9,  Figure  3.10,  shows  the  results  from  a  constant  volume,  isothermal 
system.  Relative  to  run  3,  the  initial  02  concentration  was  reduced  and  the 
initial  temperature  was  Increased.  In  Figure  3.11,  the  profiles  of  boron 
containing  species  for  4  related  cases  are  presented.  In  part  (a)  of  the 
figure,  the  initial  CO  mole  fraction  was  converted  to  H2  (i.e.,  the  system  is 
carbon  free);  in  part  (b),  the  initial  H2  mole  fraction  was  converted  to  CO; 
in  part  (c),  the  initial  02  mole  fraction  was  reduced  from  0.07  to  0.005,  by 
volume;  and  in  part  (d),  the  initial  H2  and  CO  mole  fraction  were  converted  to 
N2. 

As  a  final  set  of  examples  (Runs  11  and  12),  the  series  of  runs  9  and  10 
were  repeated,  but  with  the  initial  HBO  replaced  by  BO  +  1/2  H2 .  The  results 
from  this  series  of  runs  are  reported  in  Figures  3.12  and  3.13. 
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Figure  3.3.  Species  end  Temperature  Profiles  for  an  Adiabatic,  Constant 

S  Pressure,  Oxygen  Poor  System.  Initial  Reactants  Derived  From 

Fuel  Rich  (F/0  ■  4.5)  Equilibrium  Gas  Phase  Boron  Species. 
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3.3.  (Continued)  Species  and  Temperature  Profiles  for  an  Adiabatic, 
Constant  Prassura,  Oxygen  Poor  System.  Initial  Reactants 
Derived  From  Fuel  Rich  (F/0  "4.5)  Equilibrium  Gas  Phase 
Boron  Species. 
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Figure  3.11.  Boron  Speclation  Profiles  for  Four  Belated  HBO,  H2,  00,  02  end 
*2  Mixtures:  (a)  systea  is  carbon-free;  (b)  systeu  is  free  of 
®2 »  •7«tea  is  oxygen  deficient;  sad  (d)  systea  is  free  of 
H2  and  carbon. 
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Figure  3.12. 


Speelea  and  Temperature  Profile*  for  an  Isoeha 
Volume  Mixture  of  BO,  H2 ,  002,  02  and  N2 
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figure  3*13*  Boron  Speciation  Profiles  for  Four  Related  BO,  H2,  CO,  02 
and  M2  Mixtures:  (a)  system  is  carbon-free;  (b)  systee  is 
free  of  H2;  (c)  system  is  oxygen  deficient;  and  (d)  system 
is  free  of  H2  and  carbon. 
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Examination  of  the  concentration  and  temperature  profiles  suggest  that 
for  stoichiometric  and  oxygen  excess  systems  (l.e.,  systems  with  reaction 
products  dominated  by  HB02  and  B203),  the  reaction  proceeds  in  a  sequential, 
partially  overlapping  process  with  the  following  order  of  events  (events  here 
referring  to  the  formation  of  major  intermediates  and  final  products), 


ABC 


The  amount  of  HB02  vs  B203  in  the  products  is  determined  by  the  final 
temperature  and  the  amount  of  water  vapor  (or  more  generally,  hydrogen 
containing  species)  in  the  system,  which  is  governed  to  a  large  extent  by  the 
equilibrium  relation 

HOBO  +  HOBO  -  B203  +  ^0 

The  equilibrium  constant  for  this  reaction  is  shown  in  Figure  3.14,  and  as  can 
be  seen  the  reactants  are  generally  favored  with  the  equilibrium  shifting 
toward  the  products  with  increasing  temperature.  Over  the  temperature  range 
1800  -  2400  K,  steps  A  and  B  are  observed  to  be  fast  and  are  generally 
completed  within  a  few  micro-seconds.  The  total  energy  release  for  these  lean 
systems  occurs  over  a  lengthy  period  of  time  starting  at  a  few  micro-seconds 
and  lasting  a  couple  of  milliseconds.  Step  B  accounts  for  approximately 
33-55Z  of  the  total  energy  release  from  the  boron  containing  species.  Step  C 
(and/or  D)  accounts  for  the  remaining  50-66%  of  the  energy  release. 

For  oxygen  deficient  systems,  parallel  and  competing  steps  become 
Important  as  represented  by  the  following  network  of  events 


log  K  (c) 


Furthermore,  the  reaction  tlaes  are  substantially  longer  and  the  reaction 
products  are  dominated  by  HBO,  BO,  B202  and  HOBO  In  oxygen  deficient  systems. 

The  presence  of  GO  or  C02  can  have  a  minor  effect  on  reaction  times  and  a 
more  pronounced  affect  on  the  amount  of  B02  formed  during  the  reaction  process 
(as  seen  In  Figures  3-1 lc  and  3-1 Id  and  Figures  3-1 3c  and  3-13d). 

The  Influence  of  hydrogen  on  the  reaction  Is  shown  by  comparing  Figure 
3-13a  with  Figures  3-13d  or  13b .  In  addition  to  shortening  the  time  required 
for  equilibrium,  the  presence  of  hydrogen  plays  a  more  Important  role  In 
accelerating  the  oxidation  of  Intermediates  (i.e.,  BO  and  B02).  This  Is 
particularly  Important  with  regard  to  the  energy  release  of  the  system.  Note 
that  in  comparing  Figures  3-13a  with  3-1 3b,  the  energy  release  duration  (as 
demoted  by  the  time  to  the  peak  In  the  HB02  concentration  profile  (Figure  13a) 
or  the  time  to  the  peak  In  the  B20j  concentration  profile  (Figure  13b))  Is 
about  100  times  shorter  with  hydrogen  present. 

Comparison  of  all  the  system  with  HBO  present  Initially  with  those  with 
10  present  as  a  reactant  shows  an  Important  characteristic  difference  among 
these  two  types  of  systems.  HBO  systems  have  Induction  periods  whereas  BO 
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and  Kt'.t')  •  1.  lit)  la  da  1  a  V  laaahdaa  Mtrit  udth  iliaai 
Wf/lcj .  The  initial  aaaiUlM  (t|/l«j)|  la  da  aaca  aactar.  lailaaa 


1. 


Via  an— itivlty  aq— clone  art  seised  with 
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TIm  Genoa's  faactlaa  nutria  haa  aa  lapartaat  physical  interpretation. 

The  alaaaaaa  a#  the  Oaaaa'a  faactlaa  aatrla  nay  ha  identified  aa  the  gradients 
Ojj  •  lct(t)/»cJ(t').  t  >  t*.  These  taaaa  — y  ha  directly  Interpreted 
aa  "aaawry  feactla— “  la  that  each  tara  praaaata  tha  sensitivity  of  species  1 
at  tlaa  t  mi th  respect  ta  variations  of  species  j  at  prior  tlaa  t ' .  The 
aaaalt  la  a  aery  detailed  prehe  far  tha  aparatlaaal  aachaaiatlc  pathways  of  a 
aaapiaa  cheat  cal  aye tea.  When  t '  equals  aero,  these  gradients  are  equivalent 
ta  tha  llaaar  aaaaltivity  gradients  uhere  the  perturbed  parameter  la  the 
laltlal  age dee  eaacaat rat leas  or  initial  taaperature  of  the  ayataa. 

Fur tha raare ,  thaaa  grad lasts  are  aeaeurahle  la  the  laboratory. 

f aaapiaa  af  thaaa  laltlal  condition  aaaaltivity  gradients  are  shown  In 
Figures  3.15  threu0i  3. It.  The  eaaaplas  era  taken  froa  tha  two  kinetic 
eye  fa—  ill— traced  in  Figures  3.10  and  3.12.  Figures  3.15  and  3.17  show 
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Figure  3.13.  Sensitivity  of  S02  Concentration  Versus  TIm  to  Variations  In 
ths  Concantration  of  Several  React ants  in  an  Isothermal, 
Constant  Voltma  Mixture  of  H10,  B2,  C02,  02  and  d2 
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I  Figure  3.17.  Sensitivity  of  H  Concentration  Versus  Tine  to  Variations  in 

tha  Concentration  of  Several  Reactants  in  an  Isothermal, 
Constant  Vo lone  Mixture  of  HBO,  H2,  C02  02  and  N2 
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the  gradient*  of  9(B02  (t)J/3cj(0)  and  }[I(t)]/}Cj(0),  respectively,  for 
the  HBO  system  of  Figure  HO.  Figures  3.16  and  3.18  report  the  sane 
gradients  except  for  the  BO  system  of  Figure  3.12.  The  cbealeal  nsess  in  the 
upper  left  hand  corner  of  each  figure  denote  the  gradient  plotted.  For 
exaaqple,  the  notation  ”B02,  B202"  of  Figure  3.15(a)  represents  the  gradient 
d[B02(t)]/3[B202(0)].  The  physical  interpretation  of  this  gradient  is  slaply 
that  the  nonlnal  value  of  the  B02  concentration  is  Increased  at  a  tine  of  1  us 
and  then  decreased  slightly  thereafter  (relative  to  its  original  cone,  profile 
shown  in  Figure  3.10)  when  the  B202  concentration  is  perturbed  at  t**0.  It 
also  nay  be  Interpreted  as  a  shift  in  the  B02  concentration  profile  to  shorter 
tines,  and  thus,  an  acceleration  In  the  overall  reaction. 

Without  looking  into  the  details  of  these  sets  of  figures,  they  reveal  an 
inportant  characteristic  feature  about  the  two  systens  under  consideration. 
Note  in  particular  the  slnllarltles  in  many  of  the  gradient  profiles  of  the 
HBO  system.  This  similarity  Indicates  a  high  degree  of  dynamic  cotpling  anong 
the  various  species.  Furthermore,  the  magnitudes  of  the  coefficients  are  much 
larger  for  the  HBO  system,  and  hence,  this  system  Is  much  more  responsive  to 
perturbations.  The  high  degree  of  coupling  Is  only  in  the  vicinity  of  HBO 
consumption  and  BO  and  B02  production.  At  times  thereafter,  the  two  chemical 
systems  respond  similarly  to  Initial  condition  perturbations.  This  latter 
fact  may  be  observed  In  Figures  3.15(p)  and  3.16(m)  or  Figures  3.15(q)  and 
3.16(1)  by  comparing  profiles  at  times  greater  than  10~5  s.  The  practical 
importance  of  these  results  are  the  following:  first,  for  reaction  times 
greater  than  approximately  lus  (l.e.,  the  period  of  energy  release),  both  HBO 
and  BO  systens  respond  similarly,  and  hence,  accurate  Information  In  the 
Initial  conditions  (l.e.,  which  species  are  coming  off  the  metallic  boron 
surface)  Is  not  critical  to  the  present  study  of  the  gas  phase  chemistry.  It 
is  inportant,  however,  to  realise  that  this  statement  pertains  only  to 
stoichiometric  and  oxygen  excess  systens.  Secondly,  the  significant 
differences  In  the  responses  of  these  two  systens  (as  Indicated  by  the  Green's 
function  coefficients)  may  be  used  to  differentiate  if  HBO  or  BO  is  the 
dostfaate  species  being  produced  at  the  surface.  For  exanple,  note  that  C02 
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•iMtiiai  taa 4m  U  laUMt  the  BO  ajrtt«  (or  to  teruic  tka  B02 
mutmlw),  ■>»«—  00]  aiiltlM  to  tha  BO  tjr«t«  accalantM  it  (aid 
lacraiaaa  the  BO]  aaacaatratlaa). 

3*5.3  llaaar  FaraMtrk  Sensitivities 

Paraaatrle  sensitivities  of  fcbo  fora  dlnC2/dlaPj  bora  boon 
Investigated  for  aavaral  klooelc  system.  Scat  raaulta  fro*  thro*  of  tha 
ayataM  ora  discussed  bora.  Tbooo  lncludo  tbo  02  excess  HBO  ayataa  of  Plguro 
3*10,  tbo  02  tscaai  BO  ayataa  of  Pigurs  3.12,  and  tbo  02  doflciont  HBO  ayataa 
of  Plgura  3.11c. 

Although  on*  aay  apply  and  intarprat  tha  abora  par an* trie  derivative*  la 
aavaral  my a,  tha  noat  alanantary  Intarpratatlon  la  uaad  bora.  In  particular, 
lta  gradient  lloC]/ilnPj  will  be  Interpreted  aa  the  relative  raaponae  of 
apeciea  C2  at  tin*  t  to  a  variation  in  the  kinetic  rate  constant  par ana ter 

pj* 

Pigurea  3.19  and  3.20  show  the  gradients,  din  (B02)/31nPj,  for  the 
oxygen  excess,  HBO  and  BO  systans,  respectively.  The  sequence  of  Pigurea  (a) 
through  (p)  Includes  the  first  16  reaction  rate  paroneters  of  greatest 
laq>ortance  on  the  B02  concentration  profile.  The  figures  are  arranged  In  a 
decreasing  order  of  Importance  as  determined  by  the  auxlnun  gradient 
(absolute)  value  over  the  tine  Interval  studied.  Furthermore,  the  gradients 
plotted  are  noraallzed  to  their  maximum  values,  i.e.. 


LINEAR  SENSITIVITY  GRADIENT  LINEAR  SENSITIVITY  GRADIENT 


Figure  3.19.  Sensitivity  of  BO2  Concentration  Varans  tins  to  Variations  in 
Kinetic  Kata  Constanta  for  an  Isothermal,  Constant 
Volume  Mixture  of  HBO,  83,  CO2 ,  O2  and  O2 
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Figure  3.19*  (Continued)  Sensitivity  of  B02  Concentration  Varaua  tlaa  to 
Verletlona  In  Kinetic  Rate  Constanta  for  an  Isotberaal, 
Constant  Voluaa  Mixture  of  HBO,  H2,  C02,  02  and  02 
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Figure  3*19. 


(Continued)  Sensitivity  of  B02  Concentration  Versus  tins  to 
Variations  in  Kinetic  Rate  Constants  for  an  Isotheraal, 


Constant  Vol 


Mixture  of  HBO,  H2, 
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C02 ,  02  and  02 


LINEAR  SENSITIVITY  GRADIENT  LINEAR  SENSITIVITY 
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FI fur*  3.20.  Sensitivity  of  B02  Concentration  Versus  Time  to  Variations  In 
Kinetic  Rate  Constants  for  an  Isothermal,  Constant  Volume 
Mixture  of  BO,  H2,  C02 ,  02  and  N2 
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Figure  3.20*  (Continued)  Sensitivity  of  BO2  Concentration  Versus  Tiiae  to 
Aviations  in  Kinetic  Rate  Constants  for  an  Isothernal, 
Constant  Voluae  Mixture  of  BO,  H2 ,  CO2,  O2  and  N2 
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Figure  3.20*  (Continued)  Sensitivity  of  BO2  Concentration  Versus  Tine  to 
Variations  in  Kinetic  Rate  Constants  for  an  Isothermal, 
Constant  Volume  Mixture  of  BO,  Hj ,  CO2,  O2  and  N2 
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LINEAR  SENSITIVITY  GRADIENT 


Figure  3.20.  (Continued)  Sensitivity  of  BO2  Concentretlon  Versus  Tine  to 
Vsrletlons  In  Kinetic  Rets  Constants  for  an  Isothermal, 
Constant  Volume  Mixture  of  BO,  H2 ,  CO2,  O2  and  N2 
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Omliirln  the  MO  ijntas,  the  reactloes  seat  laportaat  to  the 
iectlao  of  Mj  are  the  breaching  roactioe 

V 

■0  ♦  Oj  ♦  *>2  ♦  0. 

(11) 

1  *• 

dleeociatioa  reaction 

- 

BBO  +  M-tB  +  BO  +  M, 

(12) 

■ 

md 

the  three  abetraetioa  reactions 

1  ■ 

HBO  ♦  0  ♦  BO  ♦  OB  (broaching) 

(13) 

■ 

■BO  ♦  OB  ♦  BO  «■ BjO 

(14) 

i 

BIO  ♦  B  ♦  10  ♦  B2 . 

(15) 

The 

reactions  laportant  to  B02  destruction  are. 

1 

BO  +  HjO  +  HB02  ♦  OH 

(16) 

| 

I02  «*■  K>2  *  +  0 

(17) 

B02  +  OH  ♦  HB02  +  o 

(18) 

8  and 

B02  +  Hj  ♦  HB02  +  H  (not  shown). 

(19) 

The 

m 

CO  oxidation  reaction 

1 

B02  +  CO  ♦  BO  +  C02 

(110) 

|  also  has  significant  influence  on  the  B02  concentration  profile. 

1  B°2 

Froe  the  C0/H2/02  reaction  subaechanlsu,  the  reactions  lnportant 
concentration  profile  are  the  initiation  reaction 

to  the 

1 

H2  +  02  ♦  H02  +  H 

(111) 

■ 

and 

the  branching  reactions 

1 

H  +  O2  OH  +  0 

(112) 

■ 

O  +  l^O+OH  +  OH  (not  shown) 

(113) 

8 

and 

H02  +  H  •*>  0  +  OH  +  M. 

(114) 

i 

0 
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la  cut  ram  t,  the  reeulta  of  Figure  3.20,  l.o.,  the  >0  system,  indicate 
that  the  aoat  acuities  rooctlou  to  Mj  product loo  la  tl.  However,  the  seat 
fuaetlcaa  (BB-tt)  era  important  to  IOj  couumptlon.  Since  BO  reacts  directly 
with  ao lecelar  o ay gee,  the  Initiation  reaction,  ill,  la  not  important,  but  the 
■|/0>2  breaching  rooctlou,  R12  and  113,  roaain  important . 

The  significance  of  backward  reactions  having  gradients  of  the  seat 
magnitude,  but  of  opposite  sign  as  their  corresponding  forward  reactions,  as 
shown  in  several  parts  of  Figures  3.19  and  3.20,  are  an  indication  of  the 
importance  of  these  reactions  equilibrium  constants  on  the  concentration 
profile  of  interest.  For  example.  Figures  3.20(b)  and  3.20(c)  show  that  for 
times  greater  than  approxieataly  100  ps,  the  equilibrium  constant  of 
B02  +  HjO  -  HB02  +  OH  is  more  important  on  the  B02  concentration  profile  than 
either  the  forward  or  backward  ruction  rate  constants. 

Sxamples  of  gradients  for  the  H-aton  concentration  profiles  for  the  three 
kinetic  systems  of  Figure  3.10,  3.12  and  3.11(c)  are  shown  in  Figure  3.21, 
3.22,  and  3.23,  rupectlvely.  Comparison  of  these  gradients  with  those  of 
Figures  3.19  and  3.20  illustrates  that  the  ranking  of  reactions  can  change 
with  each  species  of  Interest.  However,  examination  of  the  gradients  for  all 
species  reveal  a  general  grouping  of  important  reactions.  The  results  of  this 
grouping  are  reported  in  the  next  section. 

k  detailed  examination  of  the  linear  sensitivity  gradients,  as  depicted 
in  Figure  3.19  through  3.23,  reveals  an  unlimited  amount  of  information  about 
the  underlying  kinetic  processes.  For  example,  note  from  Figure  3.22((a)  that 
the  formation  of  HBO  in  the  excess  oxygen  BO  system  results  from  the 
initiation  step 

BO  +  ^  HBO  +  H  (R15) 


(R?) 


The  HBO  is  rapidly  converted  back  to  BO  mainly  through  the  propagation 
reaction 


a/2 
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Figure  3.21.  (Continued)  Sensitivity  of  H  Concentration  Versus  Time  to 
Variations  in  Kinetic  Rate  Constants  for  an  Isothensal, 
Constant  Volume  Mixture  of  HBO,  Hydrocarbon  Coabustion 
Products  and  Excess  Oxygen 
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Figure  3.21. 


(Continued)  Sensitivity  of  H  Concentration  Versus  Time  to 
Verlatlons  in  Kinetic  Rate  Constants  for  an  Isothermal, 
Constant  Volume  Mixture  of  HBO,  Hydrocarbon  Combustion 
Products  and  Excess  Oxygen 

3-73 


y; y  ;y  ;• . 


yy'.’l 


LINEAR  SENSITIVITY  GRADIENT 
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Figure  3*21.  (Continued)  Sensitivity  of  H  Concentration  Versus  Tine  to 
Variations  in  Kinetic  Rate  Constants  for  an  Isothermal, 
Constant  Voluae  Mixture  of  HBO,  Hydrocarbon  Combustion 
Products  and  Excess  Oxygen 
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LINEAR  SENSITIVITY  6RA0IENT 


LINEAR  SENSITIVITY  GRADIENT  LINEAR  SENSITIVITY 
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Figure  3*22*  (Continued)  Sensitivity  of  H  Concentration  Versus  Tins  to 
Variations  In  Kinetic  Kate  Constants  for  an  Isothermal, 
Constant  Volvass  Mixture  of  BO,  Hydrocarbon  Con  bus  ti  on 
Products  and  Excess  Oxygen 


3-76 


LINEAR  SENSITIVITY  GRAOIENT  LINEAR  SENSITIVITY  GRAOJ 
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Figure  3.22.  (Continued)  Sensitivity  of  H  Concentration  Versus  Ties 
Var*at*®**  *n  Kinetic  Kate  Constants  for  an  Isothermal 
Constant  Voluas  Mixture  of  BO,  Hydrocarbon  Combustion 
Products  sad  Excess  Oxygen 
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Figure  3.22.  (Continued)  Sensitivity  of  H  Concentration  Versus  Tine  to 
Verlations  In  Kinetic  Kate  Constants  for  an  Isothermal, 
Constant  Voluaa  Mixture  of  BO,  Hydrocarbon  Combustion 
Products  and  Excess  Oxygen 
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Figure  3.22.  (Continued)  Seneltlvley  of  H  Concentration  Versus  Tine  to 
Variations  in  Kinetic  Rate  Constants  for  an  Isothermal, 
Constant  Volume  Mixture  of  BO,  Hydrocarbon  Combustion 
Products  and  Excess  Oxygen 
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LINEAR  SENSITIVITY  6RA0IENT  LINEAR  SENSITIVITY  GRADIENT 
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Figure  3*22*  (Continued)  Sensitivity  of  H  Concentration  Versus  Tiae  to 
Verlatlons  In  Kinetic  Rate  Constants  for  an  Isothermal, 
Constant  Voluae  Mixture  of  BO,  Hydrocarbon  Combustion 
Products  and  Excess  Oxygen 
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LINEAR  SENSITIVITY  GRADIENT  LINEAR  SENSITIVITY  6RADIENT 


LINEAR  SENSITIVITY  GRADIENT 


Figure  3*23.  (Continued)  Sensitivity  of  H  Concentration  Versus  Tine  to 

Variations  in  Kinetic  Rate  Constants  for  an  Oxygen  Deficient 
Mixture  of  HBO  and  Hydrocarbon  Combustion  Products 
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LINEAR  SENSITIVITY  GRADIENT 
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Figure  3.23 
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(Continued)  Sensitivity  of  H  Concentration  Versus  Tiae  to 
Variations  in  Kinetic  Rate  Constants  for  an  Oxygen  Deficient 
Mixture  of  HBO  and  Hydrocarbon  Coabuation  Products 
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Figure  3.23. 


(Continued)  Sensitivity  of  H  Concentration  Versus  Tine  to 
Variations  in  Kinetic  Rate  Constants  for  an  Oxygen  Deficient 
Mixture  of  HBO  and  Hydrocarbon  Combustion  Products 
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Or  not*  that  la  tha  02  dafldant  ifitn,  tha  raaction  BO  +  02  ♦  B02  +  0  is 
again  an  laportant  branching  reaction. 

3*6  Critical  Inaction  Pathways 

The  critical  raaction  pathways  vara  determined  by  examining  all  the 
linear  sensitivity  gradients  Bln  C^/ainPj .  The  most  sensitive  boron 
reactions  mere  selected  as  those  which  had  gradient  values  within  2  orders  of 
magnitude  of  the  maximum  gradient  value  for  a  particular  species.  The  results 
of  this  ranking  for  excess  oxygen  systems  is  tabulated  in  Table  3-8.  For 
oxygen  deficient  systems,  the  B202  reactions  listed  in  Table  3-9  were  found 
important  in  addition  to  those  reactions  listed  in  Table  3-8. 

3.7  Extension  of  Sensitivity  Analysis  Techniques 

Currently  available  sensitivity  snalysls  techniques  have  already  been  put 
to  significant  use  In  the  first-phase  of  our  boron  combustion  study  (see 
Section  3.5  and  3.6  for  details).  The  present  studies  have  focused  on 
gas-phase  chemistry  alone  and  we  have  therefore  been  able  to  directly  utilise 
existing  computer  coding  and  experience  gained  from  sensitivity  and  modelling 
studies  of  an  allied  nature.  A  complete  model  for  particulate  boron  assisted 
hydrocarbon  combustion,  however,  must  Incorporate  heterogeneoue  and  transport 
processes.  In  addition,  the  complexities  of  the  physical  and  chemical  events 
makes  it  highly  desirable  to  search  for  any  simplifying  clues  hidden  in  the 
kinetic  behavior.  This  latter  issue  is  directly  connected  with  the 
possibility  of  scaling  and  self  similarity  behavior  existing  amongst  the 
strongly  coupled  kinetic  species  in  the  boron  particle  combustion  system. 
Therefore,  in  the  first  phase  of  this  project  we  have  examined  two  topics:  a) 
the  need  for  modifications  of  existing  sensitivity  analysis  techniques  for 
particle-gas  phase  transport  and  surface  kinetics;  b)  the  development, 
analysis  and  physical  understanding  of  scaling  and  self-similarity  behavior  in 
kinetic  systems.  The  current  status  of  these  two  items  is  summarised  below. 
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Table  3-8  -  Critical  Elementary  Reactions  for  Oxygen  Rich  Combustion 
Homogeneous 


HBO/BO  MRACTTOMS 


HBO  ♦  OH  -  BO  +  HjO  (important  when  HBO 

HBO  +  0  -  BO  +  OE  is  initial  reactant) 

HBO  +  0  B02  +  H 

HBO  +  H  -  BO  +  H2 
HBO  +  M+  BO+H  +  M 


BO  ♦  OH  -  B02  +  H  (important  when  BO 

is  initial  rectant) 

BOo  REACTIONS 

BO  +  O2  -  B02  +  0 
B02  +  00  -  BO  +  C02 


HBO2  REACTIONS 

B 02  +  H  +  M  -  HOBO  «•  M 
B02  +  H20  -  HOBO  +  M 
B©2  +  OH  -  HOBO  +  0 
B02  +  H2  -  HOBO  +  H 

HBO  +  OH  -  HOBO  +  H  (important  when  HBO 

is  initial  reactant) 


B-,0?  REACTIONS 


B02  +  B02  —  B2  O3  +  0 


Table  3-9  -  Critical  Reactions  for  Fuel  Rich  Homogeneous  Combustion 


I 

I 


B?0?  REACTIONS 

BO  +  HBO  •  B202  +  H 
BO  +  BO  +  M  -  B202  +  M 
BO  +  HOBO  -  B202  +  OH 
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3.7.1  Sensitivity  Analysis  Techniques  as  Applied  to  Boron 

Particle-Gas  Phase  Transport  and  Heterogeneous  Surface  Reactions 

Many  of  the  general  sensitivity  analysis  techniques  needed  to  treat  the 
various  aspects  of  boron  particle  combustion  are  already  available  from  prior 
applications  to  other  kinetics  problems.  The  one  exception  to  this  situation 
may  be  encountered  when  condensation  phenomena  are  ultimately  included.  Since 
modeling  of  the  product  condensation  phase  of  boron  combustion  will  only  occur 
at  an  advanced  stage  of  the  research  it  will  not  be  discussed  here.  Instead 
we  will  focus  on  the  initial  heterogeneous  burning  phase.  Currently,  the 
basic  two  problems  are  1)  transport  of  chemical  species  and  energy  to/from  the 
particle  surface,  and  2)  reaction  and  transport  on  the  particle  surface. 
Regarding  both  types  of  problems  It  is  important  to  note  that  other  steady  and 
unsteady  reaction  diffusion  problems  Including  sensitivity  analysis  have 
already  been  treated.  The  primary  distinction  is  that  these  latter  studies 
have  thus  far  been  restricted  to  one  spatial  dimension,  although  there  Is  no 
Inherent  reason  (except  for  the  level  of  computational  effort)  constraining 
them  in  this  way.  One  immediate  application  even  in  one  dimension  arises  by 
treating  the  boron  particle  burn-up  problem  as  purely  radial.  Relaxation  of 
this  Idealization  would  require  the  introduction  of  an  angular  variable.  The 
basic  mathematical  formulation  is  already  in  hand,  but  careful  consideration 
to  computer  coding  and  efficiency  would  have  to  be  made.  In  the  case  of 
adsorption,  desorption,  diffusion  and  reaction  occurring  on  the  boron  particle 
surface,  the  physical  and  mathematical  formulation  is  quite  similar  to  the 
latter  gas-phase  processes.  The  surface  problem  is  further  compounded  by  the 
fact  that  the  surface  may  be  heterogeneous  with  defects.  An  important  caveat 
is  that  the  proper  parametric  values  for  virtually  all  of  the  surface 
processes  are  far  more  speculative  than  in  the  gas-phase.  This  situation 
calls  for  special  care  when  interpreting  the  results,  but  also  provides  an 
even  more  serious  role  for  sensitivity  analysis  as  a  means  to  guide  further 
modeling  and  experimental  efforts. 
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The  gas-phase  and  surface  processes  are  intrinsically  coupled  although  a 
first  cut  at  these  probleas  would  treat  then  separately  involving  source 
^fluxes  to/fron  each  phase.  For  exanple ,  the  current  gas-phase 
aodellng/sensltlvlty  effort  could  easily  be  extended  to  include  flux  Inputs 
and  outputs  of  species  and  energy  to/fron  the  boron  particle.  Similarly,  the 
processes  on  the  particle  surface  could  also  be  treated  with  gas-phase  fluxes 
as  source/slnk  terns.  The  full  coabination  of  both  coupled  system  night  still 
be  fruitfully  treated  by  splitting  techniques  whereby  each  aspect  of  the 
problea  is  alternatively  solved  in  a  bootstrap  type  process.  All  of  the 
particulars  of  these  issues  will  be  explored  and  lnplemented  in  the  lnsulng 
research. 

3.7.2  Scaling  and  Self  Similarity  Relations  in  Kinetics 

A  kinetic  system  with  N  chemical  species  and  M  parameters  may  exhibit  a 
variety  of  complex  behavior  regardless  of  whether  it  is  steady  or  unsteady. 

In  addition,  the  nuaber  of  basic  questions  or  interrelationships  between  the 
various  variables  will  give  rise  to  N(N+M)  sensitivity  coefficients  as  a 
measure  of  the  systea  coupling.  Each  of  these  coefficients  in  turn  will  be  a 
function  of  the  coordinates  and/or  time.  As  we  have  already  seen  in  the 
present  first  phase  of  the  project  the  analysis  of  these  coefficients  can 
provide  a  critical  guide  to  physical  understanding  and  intelligent  planning  of 
subsequent  experiments.  From  a  number  of  previous  kinetic  studies  an 
intriguing  prospect  of  significantly  simplifying  this  problem  has  arisen.  In 
particular,  evidence  exists  for  the  presence  of  scaling  and  self  similarity 
relations  amongst  all  the  systea  sensitivity  coefficients  under  appropriate 
conditions. 

The  scaling  relations  take  on  the  following  form  for  a  purely  temporal 
problea 
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nhan  8j  is  a  system  parameter  such  as  rats  constant  and  Cg  li  a  chanical 
spaelaa  concantration.  This  scaling  relation  shows  that  knowledge  of  the 
sensitivities  to  one  cheaical  species  nap  be  used  to  deduce  others  by  siaply 
scaling  it  with  a  ratio  of  the  species  fluxes.  Sven  wore  remarkable  are  the 
self  similarity  conditions  which  haws  the  fore 


The  functions  Xn(t)  are  characteristic  of  the  chemical  mechanise  and  can 
most  easily  be  identified  from  the  tenporal  behavior  of  the  sensitivity 
coefficients.  In  contrast,  the  coefficients  aj  ere  constant  in  time  but 
otherwise  complicated  functions  of  all  the  systems  parameters.  The  self 
similarity  relation  is  consistent  with  the  scaling  relation  but  goes  further 
to  imply  that  knowledge  of  one  vector  of  constants  a  and  a  single  function 
XB(t)  will  yield  all  of  the  system  sensitivity  coefficients. 

The  relations  above  are  only  approximate  and  a  critical  issue  is  to 
establish  their  realms  of  validity.  The  evidence  available  thus  far  implies 
that  these  relations  are  nost  valid  when  the  system  is  characterised  by  a 
strong  nonlinear  coupling,  such  as  through  the  temperature.  Some  of  the 
current  gas-phase  boron  combustion  sensitivity  analysis  calculations  show 
evidence  of  these  relations  and  the  natter  will  be  explored  further  In  the 
second  phase  of  the  research.  The  primary  motivation  for  this  effort  is  to 
fully  utilise  any  systematic  behavior  found  in  the  kinetic  system.  An 
important  spin-off  of  these  results  concerns  their  potential  generality  since 
they  should  be  applicable  to  a  variety  of  other  siallar  problems. 
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Aa  la  th«  prutot  first  phttt  study.  It  is  antldpstsd  that  sensitivity 
techniques  will  continue  to  provide  e  powerful  aeens  to  analyse  the  physical 
■ode Is.  In  eddltion  as  the  ezperlnental  aspects  of  this  project  ult last sly 
develop  sensitivity  analysis  will  be  used  In  a  recursive  way  to  assess  the 
theoretical  and  physical  implications  the  now  experimental  information. 
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Iht  aodtl  results  sad  sensitivity  snalysls  prsssntsd  in  Ssctlon  3  fulfill 
the  major  goal  of  our  first  phase  effort;  the  basic  dependencies  of 
homogeneous  boron  suboxide  end  boron  oxyhydrlda  oxidation  are  clearly 
elucidated  and  critical  reaction  paths  (and  their  associated  rate  constants) 
and  key  laterconversion  equilibria  are  highlighted  for  further  refinement. 

Our  secondary  goal  of  identifying  issues  to  be  addresaed  in  extending  the 
model  to  the  heterogeneous  ignition  and  product  condensation  phases  was  also 
addressed. 

Among  the  major  research  conclusions  are: 

1)  The  presence  of  hydrogen  containing  species  (H20,  OH,  H,  H2)  has  a 
significant  impact  on  both  the  identity  of  species  evolved  from 
reacting  boron  particles  and  on  their  subsequent  oxidation  rata,  as 
well  as  the  rate  of  heat  release. 

2)  The  initial  Identity  of  the  boron  oxide  and  oxyhydride  species 
evolved  from  the  particulate  boron  influences  the  oxidation 
chemistry  for  a  very  short  time  period  (typically  on  the  order  of  a 
alcrosecond);  subsequent  oxidation  chemistry  at  reasonable  0/F 
ratios  is  only  mildly  dependent  on  the  identity  of  Initially 
vaporised  species  including  B,  BO,  HBO  and  B202.  This  fact 
validates  our  decision  to  nodal  the  homogeneous  combustion  phase 
chemistry  without  being  able  to  fully  specify  the  kinetic 
distribution  of  initially  volatilised  boron  species  from  Ignited 
particles. 
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3)  As  nvultd  by  our  sensitivity  analysis,  «  relatively  small  set  of 
kinetic  and  equilibrium  processes  appear  to  dominate  the  boron 
oxidation  pathways.  As  discussed  in  Subsection  4.2,  this  allows  the 
identification  of  key  kinetic  and  thermochemical  parameters  which 
need  further  experimental  and/or  theoretical  refinement  if  a  more 
accurate  model  is  desired. 

4)  The  most  important  issues  associated  with  extending  our  boron 
combustion  model  and  concomitant  kinetic  sensitivity  analysis 
techniques  to  the  initial  particle  ignition  and  heterogeneous  burn 
phase  of  boron  combustion  have  been  identified.  These  include 
transport  of  chemical  species  and  energy  to/ from  the  particle 
surface  and  kinetic  processes  on  the  surface.  This  initial  analysis 
suggests  that  both  processes  may  be  treated  with  existing 
sensitivity  analysis  techniques  by  invoking  two  simplifications. 
First,  the  gas  phase  and  surface  processes  would  be  treated 
separately  (as  in  the  present  work),  but  source  and  sink  terms 
corresponding  to  species  and  energy  fluxes  would  be  incorporated 
into  the  kinetic  equations.  Second,  reaction  and  transport  to  and 
from  the  particle  surface  would  be  modeled  for  one  spatial  dimension 
(e.g.,  the  radial  dimension).  This  would  provide  an  initial  model 
for  the  heterogeneous  burn,  a  more  realistic  model  for  the  gas  phase 
oxidation  kinetics,  and  build  in  a  coupling  between  the  two  phases 
using  chemical  species  fluxes.  On  the  basis  of  this  work,  it  would 
be  possible  to  evaluate  the  validity  of  scaling  and  self-simllarly 
relations. 

5)  A  better  defined  (but  possibly  incomplete)  model  of  homogeneous 
boron  assisted  hydrocarbon  combustion  now  exists  and  can  be  used  to 
study  a  number  of  energetic  combustion  and  propulsion  systems  of 
interest  to  the  Department  of  Defense  coemunlty.  Several  of  these 
are  discussed  in  Subsection  4.3. 
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4.2  Raco— a ndationa  for  Further  Research 


The  work  described  In  this  report  hss  significantly  increased  our  current 
understanding  of  the  hoaogeneous  gas  phase  kinetics  for  boron-assisted 
hydrocarbon  coabustlon.  However,  achieving  a  reliable  model  for  a  boron  fueled 
coebustor  will  require  further  research  on  several  fronts.  This  includes  the 
development  of  models  for  the  Initial  heterogeneous  burn  and  final 
condensation  stages  of  boron-assisted  combustion,  in  conjunction  with 
experimental  work  which  can  provide  accurate  values  for  key  thermocheaical  and 
kinetic  parameters,  as  well  as  validation  of  the  various  combustion  models. 

In  subsections  4.2.1  and  4.2.2  below,  we  briefly  discuss  our 
recommendations  for  further  experimental  and  modeling  work.  It  is  Imperative, 
however,  to  note  that  real  progress  in  this  field  requires  that  the  modeling 
and  experimental  work  be  closely  Interrelated.  In  particular,  modeling  of 
either  the  initial  ignition  and  heterogenous  oxidation  or  the  final 
condensation  processes  will  rely  heavily  on  the  homogeneous  combustion  model 
presented  in  this  report.  Consequently,  any  further  development  in  the 
combustion  modeling  must  be  carefully  coordinated  with  experimental  studies 
which  can  be  used  to  validate  the  present  homogeneous  combustion  model. 
Similarly,  as  key  reaction  mechanisms  in  the  other  stages  of  boron  combustion 
are  identified,  additional  experimental  work  will  be  required  to  provide 
accurate  values  for  the  input  parameters  and  for  comparisons  with  model 
predictions. 

4.2.1  Further  Basic  Experimental  and  Theoretical  Research 

In  the  present  program,  sensitivity  analysis  was  used  to  identify  the 
critical  chemical  and  physical  parameters  which  control  the  homogeneous 
combustion  process.  The  capability  to  identify  the  most  important  gas  phase 
reactions  allows  researchers  to  preferentially  select  those  rate  parameters 
whose  magnitudes,  temperature  dependence,  and  product  branching  ratios  need  to 
be  carefully  evaluated  with  further  experimental  work.  Such  experimental 
studies  are  particularly  important  for  modeling  of  boron-assisted  fuels 
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because  of  the  large  uncertainties  In  both  the  overall  boron  combustion 
process  and  specific  Input  parameters.  Thus,  while  the  estimated  rate 
parameters  from  the  present  and  earlier  work  were  adequate  for  initial 
modeling  efforts,  they  cannot  be  expected  to  be  accurate  enough  for  the 
modeling  which  will  be  required  to  optimize  boron  combustion  systems. 

Based  on  the  results  from  the  present  modeling  of  the  homogeneous 
combustion  phase.  It  Is  clear  that  experimental  measurements  of  the  kinetic 
rate  parameters  and  product  branching  ratios  for  HOBO,  HBO,  BO,  B02,  and  B202 
reactions  (Table  3-6)  are  needed.  The  kinetic  sensitivity  analysis  described 
in  Sections  3.5  and  3.6  suggests  that  experimental  studies  should  particularly 
concentrate  on  the  reectlons  listed  in  Tables  3-8  and  3-9  since  the  gas  phase 
kinetics  has  been  shown  to  be  particularly  sensitive  to  these  elementary 
reactions. 

To  date,  relevant  laboratory  studies  have  been  limited  to  reactions  of 
atomic  B(g)  and  BO(g)  with  simple  stable  oxidants  such  as  02,  C02,  and  H20.  A 
similar  situation  exists  In  for  boron/halogen  systems  where  only  a  few 
reactions  of  BF  and  BC1  with  oxidants  such  as  02,  0,  and  N02  have  been 
studied.  While  these  past  and  on-going  studies  are  Important,  only  the  work  on 
BO,  Initially  studied  by  Fontijn4"1  and  more  recently  by  Oldenborg  and 
Baughcum,**-2  is  relevant  to  boron-assisted  hydrocarbon  oxidation,  since  only 
BO,  of  the  gaseous  species  studied  to  date,  is  expected  to  be  vaporized  from 
boron  particles  burning  In  a  hydrocarbon/oxygen  combustion  system. 
Consequently,  further  progress  in  formulating  a  reliable  boron  combustion 
model  requires  extensive  experimental  research  on  key  reactions  involving  the 
other  boron  oxides  and  oxyhydrldes  which  play  important  roles  in  the 
combustion  process. 

The  sensitivity  analysis  summarized  in  Subsection  3.5  has  also  indicated 
a  need  for  more  accurate  thermochemical  data.  Equilibrium  constants  for  key 
reactions  Involving  HB02  (g)  formation  and  conversion  to  B203  (g)  are 
particularly  Important.  Specifically,  equilibrium  data  for  two  reactions 
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B02  +  H20  -  HOBO  +  OH 


HB02  +  HB02  -  B203  +  HjO 

appear  to  be  particularly  critical.  Better  valuea  for  the  equilibrium 
constants  for  these  key  processes  would  come  from  the  experimental  kinetics 
community  (by  Independent  determination  of  forward  and  reverse  rate 
constants),  the  experimental  thermodynamics  community  through  second  and  third 
law  AH*  and  AS  determinations,  or  from  the  theoretical  molecular  structure 
community  through  better  molecular  parameters  for  HBO  (g),  B2O3  (g),  and  B02 
(g).  In  general,  reference  to  Table  3-1  Illustrates  the  crucial  need  for  more 
accurate  (hopefully  to  the  0.1  kcal/mole  range)  molecular  structure  data  over 
the  temperature  range  of  300  to  3000*11. 

4 .2 .2  Further  Modeling  Work 

Coupled  to  the  experimental  studies  described  in  4.2.1,  additional 
modeling  of  both  the  heterogeneous  burning  of  a  boron  particle  and 
condensation  of  the  gas  phase  combustion  products  Is  needed  to  develop  a 
complete  model  for  the  overall  combustion  process.  The  most  logical  extension 
of  the  present  work  Is  to  model  the  gas  phase  chemistry  in  the  presence  of  a 
clean  solid  or  liquid  boron/boron  oxide  surface.  From  a  practical  viewpoint, 
this  extended  model  will  represent  the  clean  particle  burning  phase  described 
In  Section  2.2.  From  a  more  fundamental  perspective,  the  same  model  could  be 
applied  to  analyze  simple  heterogeneous  kinetic  experiments  [see  for  example. 
Ref.  2.23]. 

This  is  a  particularly  important  aspect  of  the  practical  combustion 
problem  since  earlier  work,  as  well  as  the  results  of  the  present  study 
Indicate  the  boron  oxidation  Is  dominated  by  heterogeneous,  surface  oxidation 
processes.  This  preliminary  work  would  emphasize  the  oxidation  kinetics  of 
B(s)  by  the  various  gas  phase  species  expected  in  the  boron/hydrocarbon 
combustion  system.  Subsequent  work  should  focus  on  the  boron  particle 
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Ignition  rad  the  elementary  chemical  and  physical  processes  of  boron  oxide 
film  growth  and  vaporisation.  Particular  attention  should  be  given  to  the 
competing  diffusion  processes  of  the  various  species  in  B2O3  which  control 
oxide  growth  and  the  chealcal  transformations  affecting  the  oxide  removal  at 
temperatures  lower  than  the  vaporisation  temperature  of  B2O3  [e.g.,  B203(1)  + 
H2O  (g)J.  The  modes  of  heat  transfer  should  be  an  Important  phase  of  this 
modeling  since  it  has  been  shown  experimentally  that  particle  self-heating  is 
negligible  during  the  Ignition  process. 

finally,  the  condensation  process  which  occurs  downstream  from  the 
burning  stage  should  be  treated.  This  phase  nay  be  Initially  studied 
separately  from  the  heterogeneous  particulate  boron  oxidation  and  subsequently 
coordinated  with  the  other  submodels  to  provide  a  complete  description  of  the 
full  process. 

4.3  Belated  Areas  of  Application 

The  work  presented  in  this  report  was  specifically  motivated  by  the 
current  interest  in  boron-assisted  hydrocarbon  combustion  for  air-breathing 
propulsion  systems.  However,  metallic  boron  and  borane  based  propellants  are 
also  being  actively  considered  for  a  variety  of  related  applications. 

Although  each  potential  application  has  unique  physical  and  chemical 
characteristics,  there  will  undoubtedly  also  be  strong  similarities  in  the 
basic  combustion  chemistry.  It  is  worthwhile,  therefore,  to  briefly  review  the 
relation  between  our  modeling  program  and  ongoing  research  in  other 
laboratories . 

The  use  of  boron  enhanced  fuels  in  decoy  flares  is  one  area  of 
application  of  current  interest.  Infrared  decoy  flares  deployed  from  both 
tactical  and  strategic  aircraft  must  present  a  hot  infrared  spectral  signature 
which  will  divert  a  missile  seeker  away  from  the  aircraft  nozzle  and  exhaust 
plume.  In  order  to  be  successful,  a  decay  flare  must  achieve  very  rapid  and 
efficient  combustion  while  simultaneously  providing  the  appropriate  spectral 
signature.  Current  magnesium/ teflon  flares  have  significant  deficiencies  on 
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both  accounts.  Ongoing  work  at  Aerodyne  Research,  the  Naval  Weapons  Support 
Center  and  other  laboratories  have  identified  several  boron  enhanced  fuels 
which  may  have  both  combustion  and  exhaust  spectral  signatures  that  could 
result  in  improved,  advanced  infrared  decoys.  The  knowledge  gained  in  the 
present  AFOSR  boron  combustion  modeling  program  will  be  directly  relevant  to 
the  evaluation  of  these  new  decoy  concepts. 

Another  related  area  of  interest  is  the  addition  of  borane  or 
borane/amlne  adduct  compounds  to  liquid  rocket  propellants  in  conjunction  with 
a  variety  of  oxidizers  including  nitrogen  oxides,  oxygen,  oxygen  halides,  or 
interhalogens .  Attack  warning  issues  have  stimulated  a  great  deal  of 
Department  of  Defense  sponsored  research  effort  focused  on  determining  the 
exhaust  plume  afterburning  chemical  kinetics  and  spectral  properties  of  boron 
containing  exhaust  species  from  such  "advanced”  liquid  rocket  propellants. 

dorane-based  propellants  may  also  have  significant  applications  as  burn 
rate  modifiers  for  NC  and  NC/HMX  propellants  and  ignitors  in  large  and  medium 
caliber  guns  and  telescoped  ammunition,  and  as  propellants  for  advanced  gun 
propulsion  concepts  such  as  consumable  sabot,  FILMBAL  and  the  traveling  charge 
concept.4-3  Conventional  chemical  gun  propulsion  is  intrinsically  inefficient 
for  obtaining  muzzle  velocities  greater  than  2  km/sec.  In  a  traveling  charge 
configuration,  a  projectile  is  accelerated  by  a  combination  of  gas  pressure 
and  Impulse  produced  by  a  very  fast  burning  propellant  attached  to  and  moving 
with  the  projectile.  Initial  calculations  Indicate  that  burn  rates  between 
100  m/s  and  1000  m/s  are  required  to  obtain  the  desired  muzzle  velocities. 

Much  of  the  recent  work  on  very  high  burn  rate  (VHBR)  propellant 
formulations  has  focused  on  a  series  of  borane  salts  (B1 0B1 0-2,  B12H12-2) 
known  generlcally  by  the  trade  name  of  Hivelite.  Although  these  salts  have 
been  shown  to  be  very  effective  bum  rate  promoters,  there  is  very  little 
understanding  of  the  basic  mechanisms  through  which  they  act.  It  is  generally 
accepted  that  VHBR  materials  bum  by  a  convective  mechanism  in  which  the 
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combustion  zone  propagates  by  the  Infiltration  of  hot  combustion  products 
through  pores  in  in  the  propellant  saaple.  This  is  followed  by  a  transition  to 
a  very  rapid  combustion  which  takes  place  throughout  the  whole  saaple. 

In  borane-assisted  combustion,  the  initial  ignition  and  volatilization 
combustion  phase  associated  with  metallic  boron  propellants  is  replaced  by  a 
simpler  and  much  faster  pyrolysis  of  the  borane.  However,  the  chemical 
kinetics  of  the  subsequent  combustion  phases  should  closely  parallel  the 
reaction  kinetics  discussed  in  Section  3.5.  Specifically,  as  the  borane  ion 
breaks  up  and  undergoes  some  initial  oxidation,  the  resulting  combustion  gases 
will  probably  contain  many  of  the  same  oxides  and  oxyhydrldes  which 
characterize  the  homogeneous  phase  of  metallic  boron  combustion.  Consequently, 
the  homogeneous  boron  combustion  model  described  in  Section  3.5  should  be  a 
very  useful  in  research  on  the  mechanisms  by  which  the  gaseous  combustion 
products  Influence  the  burn  rate.  The  first  step  in  generalizing  the  reaction 
list  in  Table  3-6  to  borane-assisted  propellants  would  be  to  Include  oxidation 
reactions  for  BH  fragments  and  substitution  reactions  for  the  highly  reactive 
5-coordinated  boron  atoms  in  the  deca-borane  salt.  Some  evidence  for  the 
similarity  between  borane  and  metallic  boron  assisted  combustion  can  be  seen 
in  Table  4-1 •  The  table  gives  combustion  product  concentrations  for  a 
representative  sample  of  a  VHBR-borane. 

Table  4-1  -  Combustion  Products  of  VHBR-Borane  Propellant  and  a 
Conventional  High  Impetus  Propellant 


Species 

Concentration 

HB02(g) 

0.186 

B203(g) 

0.0072 

HBO(g) 

0.092 

B203(1) 

1.88 

It  anticipated  that  as  the  present  modeling  efforts  are  extended  to  the 
heterogeneous  phases  of  metallic  boron  combustion,  the  resulting  combustion 
models  and  mathematical  framework  developed  to  extend  sensitivity  analysis  to 
include  heterogeneous  processes  will  be  extremely  useful  in  understanding  the 
basic  mechanisms  responsible  for  the  action  of  VHBR  propellants  based  on 
boranes. 
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